Cavity cooling of internal molecular motion 
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We predict that it is possible to cool rotational, vibrational and translational degrees of freedom of 
molecules by coupling a molecular dipole transition to an optical cavity. The dynamics is numerically 
simulated for a realistic set of experimental parameters using OH molecules. The results show that 
the translational motion is cooled to few /iK and the internal state is prepared in one of the two 
ground states of the two decoupled rotational ladders in few seconds. Shorter cooling times are 
expected for molecules with larger polarizability. 
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The preparation of molecular samples at ultra-low 
temperatures offers exciting perspectives in physics and 
chemistry [ij. This goal is presently pursued by several 
groups worldwide with various approaches. Two meth- 
ods for generating ultracold molecules employ photoasso- 
ciation and Feshbach resonances, and are efficiently im- 
plemented on alkali dimers Another approach uses 
buffer gases, to which the molecules thermalize [Ij. Its 
application is limited by the physical properties of atom- 
molecule collisions at low temperatures. Optical cooling 
of molecules is an interesting alternative, but, contrary 
to atoms, its efficiency is severely limited by the mul- 
tiple scattering channels coupled by spontaneous emis- 
sion, and may only be feasible for molecules which are 
confined in external traps for very long times Ele- 
gant laser-cooling proposals, based on optical pumping 
the ro- vibrational states 0, 0] and excitation pulses tai- 
lored with optimal control theory 0, 0] , exhibit efficien- 
cies which are indeed severely limited by spontaneous de- 
cay. In [1, 0] it was argued that cooling of the molecular 
external motion could be achieved by using resonators, 
by enhancing stimulated photon emission into the cav- 
ity mode over spontaneous decay. This mechanism was 
successfully applied for cooling the motion of atoms ■ 

In this Letter we propose a method for optically cool- 
ing external as well as the internal degrees of freedom 
of molecules. The method relies on the enhancement of 
the anti-Stokes Raman transitions through the resonant 
coupling with the modes of a high-finesse resonator, as 
sketched in Fig.[TJ All relevant anti-Stokes transitions are 
driven by sequential tuning of the driving laser. At the 
end of the process the molecule is in the ro-vibrational 
ground state and the motion is cooled to the cavity 
linewidth. We demonstrate the method with ab initio 
based numerical simulations using OH radicals, of which 
cold ensembles are experimentally produced [ill . Il^ . 

We now outline the theoretical considerations. We con- 
sider a gas of molecules of mass M, prepared in the elec- 
tronic ground state X, and with dipole transitions X — s- 
E. Here, E is a. set of electronically excited states, includ- 
ing higher-lying states which may contribute significantly 
to the total polarizability. We denote the ro-vibrational 
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FIG. 1: (a) A molecular sample interacts with the cavity field 
and is driven by a laser, inducing Raman transitions cool- 
ing the internal and external degrees of freedom, (b) Comb 
of resonances at which photon emission into the cavity is en- 
hanced. The grey bars symbolize the molecular lines, which in 
OH extend over several tens nm. The laser frequency (arrow) 
is varied to sequentially address several anti-Stokes lines. 



states and their corresponding frequencies by \j, £,) and 

^5"''' (j = ^1 ^5 •■Oi while the elements of the dipole mo- 
ment d are V^^^,> ^ {A e E , \d\X , ^) . These transi- 
tions are driven by a far-off resonant laser and interact 
with an optical resonator as illustrated in Fig. [TJa). In 
absence of the resonator, spontaneous Raman scattering 
determines the relevant dynamics of molecule-photon in- 
teractions. These processes depend on the center-of-mass 
momentum p through the Doppler effect and occur at 



rate r^?;^j,(p) 
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the decay rate along the transition \ A, ^") \X, ^'), and 
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with Fj" = Z^f'Tf"^?'- Here, flL.^^i" = Sl^D^^^" ■ 
fi)/^ gives the strength of coupling to the laser with 
electric field amplitude £l, polarization e^, frequency lul, 
and wavevector kL, and ^f^>i = ^[^^ ~ w^^"* -\- ujl de- 
notes the detuning between laser and internal transition. 
In the limit where the laser is far-off resonant from the 
electric dipole transitions, electronic ground states at dif- 
ferent ro-vibrational quantum numbers are coupled via 
Raman transitions. The corresponding emission spec- 
trum is symbolized by the grey bars in Fig. [Ifb). In this 
regime, coupling an optical resonator to the molecule may 
enhance one or more scattering processes when the cor- 
responding molecular transitions are resonant with res- 
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onator modes, symbolized by the black bars in Fig. [Ijb). 
This enhancement requires that the rate r^_,^, (p), de- 
scribing scattering of a photon from the laser into the 
cavity mode and its subsequent loss from the cavity, ex- 
ceeds the corresponding spontaneous Raman scattering 
rate, F^^^, (p) ^ F^^^,(p). For a standing wave cav- 
ity, in the regime where the molecular kinetic energy ex- 
ceeds the cavity potential, and when the cavity photon 
is not reabsorbed but lost via cavity decay 0, we have 
r^_^,(p) = F^:^^,(p) + F^::^,(p), where the sign ± gives 
the direction of emission along the cavity axis and 
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^„ ((5w±kc •p/M)^ + k2 
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Here, 2k is the cavity linewidth, Suj ^ uj'^ ' — uj^, , 
lul — f^c is the frequency difference between initial and 
final (internal and cavity) state, with flc the frequencies 
of the cavity modes, and g^^,^^,, are the Fourier com- 
ponents at cavity-mode wave vector ±|kc| of the cou- 
pling strength to the empty cavity mode, ^c^'^f" (x) = 
fc(x)(eo ■ with £c and eo vacuum amplitude 

and polarization. Note that reabsorption and sponta- 
neous emission of the cavity photon, can be neglected 

while K > |5c,{'^«"^i,?— 

Enhancement of Rayleigh scattering into the cavity 
is achieved by setting the laser on resonance with one 
cavity mode, lol = ^c, see Fig. [ijb), provided that 
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:$> 1, where F is deter- 
mined by the linewidths of the excited states which sig- 
nificantly contribute to the scattering process. This sit- 
uation has been discussed in where it has been pre- 
dicted that the motion can be cavity cooled to a tem- 
perature which is in principle only limited by the cavity 
linewidth [l^, provided that the laser is set on the low- 
frequency side of the cavity resonance. Note that cooling 
of the motion in the plane orthogonal to the cavity axis is 
warranted when the laser is a standing wave field, which 
is simply found in our model by allowing for the absorp- 
tion of laser photons at wave vector — kL. In general, 
enhancement of scattering along the Raman transition 
^ — > decreasing the energy of the ro- vibrational de- 
grees of freedom, is achieved by setting the laser such 
that the corresponding anti-Stokes spectral line is reso- 
nant with one cavity mode, and requires j_,j//F«; 2> 1. 
The cooling strategy then consists of choosing a suitable 
cavity and of sequentially changing the laser frequency, 
so to maximize the resonant drive of the different anti- 
Stokes spectral lines, and thereby cooling the molecule 
to the ro-vibrational ground state. 

We simulate the cooling dynamics for OH radicals us- 
ing a rate equation based on the rates H]) and ((2]). The 
considered Raman process is detuned from the excita- 
tion energy of 32402 cm" ^ (971.4 THz) between the X^H, 
ground state and the electronically excited state A^E^ 
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FIG. 2: (a) Simulated Raman spectra for the first nine rota- 
tional states of OH, which are relevantly occupied at room 
temperature. (b) Potential energy surfaces of the X^IIi 
ground state and of the A^E+ excited state. The coherent 
Raman process is indicated by the arrows, (c) Ro-vibrational 
substructure. 



as indicated in Fig. HKb). The relevant ro-vibrational 
spectrum and the coupling of the molecule to the laser 
field and cavity modes were obtained by combining ab 
initio calculations for the vibronic degrees of freedom 
with available experimental data for the rotational con- 
stants. The level scheme is displayed in Fig. [D^b). The 
potential energy surfaces (PES) and the polarizabitilics 
were calculated with highly correlated quantum chem- 
ical methods: the electronic structure calculations 
were performed on the Multi Configurational SCF level 
using a single atom basis set (aug-cc-pVTZ). Rates ([T]) 
and ([2]) were evaluated using the polarizability tensors, 
defined as a^^^' = '^^^("'^("^C /^^i" [3 ^'^'^ 
calculated with linear response theory at the MCSCF 
level 3, 11 , l3l • III order to determine the internal level 
structure and transition strengths, the vibrational eigen- 
values and cigcnfunctions were evaluated with a relax- 
ation method using propagation in imaginary time plus 
an additional diagonalization step [Toj . The correspond- 
ing Placzck- Teller coefficients [20| were calculated for the 
transitions between the rotational sublevels. 

In order to obtain a concise picture of the cooling dy- 
namics, several assumptions were made without loss of 
generality. The molecules are prepared in the lower-lying 
component of the X^H electronic ground state. 
The hyperfine splitting is neglected as angular momen- 
tum conservation for the rotational Raman transitions 
inhibits transitions between the hyperfine sublevels. At 
T « 300 K only the first 9 rotational states of OH are rel- 
evantly occupied, while only the vibrational ground state 
is populated. The selection rule A J = 0, ±2 for the 
rotational transitions yields that the matrix elements of 
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transitions between rotational states with opposite parity 
vanish, resulting in two separate ladders for the scatter- 
ing processes with final states \X,v = 0, J = 0, 1) plj . 
We assume that a preparation step has occurred, bring- 
ing the motional temperature to below 1 K. This could 
be realized with, e.g. , helium-bufFer-gas cooling [21 , elec- 
trostatic filtering [25 



or decelerator techniques [22 
The high-finesse optical cavity has a free-spectral range 
(FSR) of 15 X 27r GHz, which can be realized with a Fabry- 
Perot-type cavity of length L = 1 cm. For simplicity we 
assume that the cavity only supports zeroth-order trans- 
verse modes. This actually underestimates the possibili- 
ties of scattering light into the cavity, since higher-order 
transverse modes can be combined in degenerate cavi- 
ties, like confocal resonators. The cavity half-linewidth 
is set to K = 75 X 27rkHz, and the coupling gc.o^o = 
2tt X 116 kHz. This is achieved with a mode volume of 
3.2 X 10^^^ m'^, assuming a mode waist of wq = 6 /im, 
and a cavity finesse F = 10^, i.e., a mirror reflectivity of 
0.999969. We also choose a laser wavelength of 532 nm, 
for which ample power is available as well as mirrors of 
the required quality. The frequency of the laser is far be- 
low that of the OH A-X ro-vibronic band. We assume to 
have single-frequency light of 10 W enhanced by a factor 
of 100 by a build-up cavity in a TEMqo mode, corre- 
sponding to a Rabi coupling JIl.o^o = 27r x 69 GHz and 
frequency lul = lu^"^^ - lu^'^^ - A with A « 27r x 407 THz. 
The latter value is sequentially varied during cooling, in 
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FIG. 3: Reduced spectrum of the relevant rotational anti- 
Stokes transitions: the lines are projected onto a single free- 
spectral range of the cavity, whose width is indicated by the 
arrows. The unfolded Raman spectrum spans 15 THz or up to 
10^ cavity modes. The frequencies of the lines are evaluated 
from quantum chemical calculations, and must be understood 
as a qualitative picture; high-resolution experimental input is 
needed to fix the absolute position with kHz accuracy. 

order to drive (quasi-)resonantly the cooling transitions. 
In combination with the broad spectrum of cavity modes, 
the laser only needs to be varied over one FSR to address 
all anti-Stokes lines. Fig. [3| displays the anti-Stokes Ra- 
man lines as a function of the frequency modulus the 
FSR. Addressing the Stokes lines (not drawn) can be 



avoided, given the small laser and cavity lincwidths. In 
a confocal cavity, e.g., all higher-order transverse cavity 
modes will be degenerate with fundamental ones. In ad- 
dition, our scheme is robust against a small number of 
coincidences between Stokes and anti-Stokes lines. 

The cooling strategy is as follows. First, the external 
degrees of freedom are cooled to the cavity linewidth, cor- 
responding to a final temperature T « 4/xK, by setting 
the Rayleigh transition quasi-resonant with one cavity 
mode. The corresponding coefficients have been evalu- 
ated numerically, giving the rate of Rayleigh scattering 
for OH into the cavity Fq^q 1 kHz, while the spon- 
taneous rate F^^^ ~ 0.5 — 2.5 Hz. We verified the effi- 
ciency of cooling by sol ving the semiclassical equation for 
the mechanical energy [25|. For these parameters, start- 
ing from T ^ IK for the external degrees of freedom, 
the cooling limit is reached in a time of the order of 1 s. 
Then, the rotational degrees of freedom are cooled by set- 
ting the laser frequency to sequentially address each anti- 
Stokes spectral line. A manually optimized sequence led 
to the result in Fig. [U where the mean rotational quan- 
tum number (J) is plotted as a function of time. The 
final value (J) « 0.5 corresponds to the final situation in 
which the two states J = and J = 1, ground states of 
each ladder, achieve maximum occupation, equal to 50%. 
The insets in Fig.[4|show that after 0.3 s their occupation 
is about 40%, while after 1.8 s it reaches 49% (leading to 
a total population of 98.8%). The cooling rate for the ro- 
tational degrees of freedom of OH is of the order of 4 Hz, 
see Fig. [H while the rate of heating due to spontaneous 
Raman scattering along the Stokes transitions is about 
0.1 Hz. In the simulation the vibrational degrees of free- 
dom are taken into account but no vibrational heating 
is observed. In a separate simulation, we checked that 
vibrational excitations are cooled with the same scheme. 

The cooling time can be improved significantly for 
molecules with higher polarizabilities a, as the cooling 
rate scales with a^. For instance, the polarizability of 
NO is approximately 4 times larger than for OH, and 
preliminary results show that it indeed cools down faster, 
although the cooling time is not reduced by a factor 16 
because more rotational states are occupied at 300 K. For 
molecules like Cs2, a is two orders of magnitude larger 
than for OH, and should yield faster cooling rates. 

For molecules with low polarizability, like OH, reduc- 
tion of the cooling time seems possible by further op- 
timizing the sequential procedure. In addition, the ef- 
ficiency could be improved by using degenerate cavity 
modes or by superradiant enhancement of light scatter- 
ing, sustained by the formation of self-organized molec- 
ular crystals [27|. 

In summary, we presented a strategy for cooling ex- 
ternal and internal degrees of freedom of a molecule. 
We simulated the cooling dynamics for OH using experi- 
mentally accessible parameter regimes, showing that this 
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FIG. 4: Mean rotational quantum number versus time during 
the cooling process. The manually optimized cooling sequence 
first empties the levels J — (2,3), thereafter higher levels 
are addressed. The small figures show the state distributions 
at time 0, 0.3 and 1.8 s. The cavity length is fine-tuned to 
address the J2^o and the J3_i transition simultaneously. 



method allows for efficient preparation in the lowest rovi- 
brational states, while the motion is cooled to the cavity 
lincwidth. For OH the cooling time is of the order of sec- 
onds, and requires thus the support of trapping technolo- 

gies which are stable over these times 0,11 [ii,!!!,!!^]. 

The cooling time of molecules with larger polarizabili- 
ties can scale down to few ms, when the polarizability is 
about 10 times larger. Applications of this technique to 
polyatomic molecules has to deal with an increasing num- 
ber of transitions to be addressed, which will slow down 
the process. A possible extension of this scheme could 
make use of excitation pulses, determined with optimal 
control techniques 01 ■ 
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